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Abstract 

Quasi-Spherical harmonics, Yf l (9, (j)) are derived and presented for 
half-odd-integer values of I and m. The form of the 4> factor is identical 
to that in the case of integer £ and m: exp (irrxft). However, the domain 
of these functions in the half-odd-integer case is 0<</><47r rather than 
the domain 0<4><2tt in the case of integer I and m (the true spherical 
harmonics). The form of the 9 factor, pj m '(#) (an associated Legendre 
function) is (as in the integer case) the factor (sin#)' m ' multiplied 
by a polynomial in cos# of degree {£— |m|) (an associated Legendre 
polynomial). A substantial difference between the spherical (integer 
£ and m) and quasi-spherical (half-odd-integer £ and m) Legendre 
functions is that the latter have an irrational factor of Vsin# whereas 
the factor of the truly spherical functions is an integer power of sm6. 
The domain of both the true and quasi spherical associated Legendre 
functions is the same: 0<9<ir. A table of the Associated Legendre 
Functions is presented for both integer and half-odd-integer values of 
£ and m, for \m\ = 0, |, 1 . . . ^, and for (^-|m|) = 0, 1, 2, 3, 4, 5. The 
table displays the similarity between the functions for integer £ and m 
(which are well known) and those for half-odd-integer £ and m (which 
have not been recognized previously). 



*Author to whom correspondence should be addressed; email: ghunter@yorku.ca 
^Permanent Address: Institute for Nuclear Energy Research, Bulgarian Academy of 
Sciences, Sofia, Bulgaria 



1 



1 Introduction 

The theory of angular momentum based upon the fundamental commutation 
relations fl], p. 93] |2|, p. 309], || pp. 107-112], produces the eigenf unctions of 
L 2 (the square of the total angular momentum) and L z (its ^-component) 
with eigenvalues of Ti 2 £(£ + 1) and km, respectively [p], §5.4 pp. 100-103]. 
This general theory leads (via the raising and lowering ladder operators) to 
the prediction of both integer and half-odd-integer values of the quantum 
numbers £ and m [[I], §5.4 pp. 100-103], p. 311], with the manifold of 
eigenvalues defined by: 



The quantum numbers, £ and m , are half-odd-integers for the spin and total 
(spin+orbital) angular momentum of fermions, and are integers for fermion 
orbital angular momentum, and for boson spin and total angular momentum. 

This general theory of angular momentum is abstract in the sense that 
the operators, L 2 and L z , and their eigenfunctions, are not functions of any 
coordinates. 

In the Schrodinger wave mechanics of the hydrogen atom the operators, 
L 2 and L z (that represent the motion of the electron around the proton) are 
expressed in terms of spherical polar coordinates: r, 8, <fi [|], p. 207], 
fl], p. 95], where they have the form: 



The eigenfunctions in this coordinate representation are called spherical har- 
monics and denoted by Y™(6, <p). This Schrodinger representation of angular 
momentum leads to the same manifold of eigenvalues as the general theory, 
except that £ and m are restricted to be integers 0, pp. 313-315]. This 
restriction arises from the argument that Y™(9, 0) should be a single- valued 
function of the coordinates |], p. 103]. 

The purpose of this paper is to derive and present the eigenfunctions of 
the coordinate-representation operators (|1|) corresponding to half-odd-integer 
values of £ and m. Their 9 factors are compared with those of the well-known 
eigenfunctions with integer £ and m in Table [I] on page ^. 



. £ = | 1 | 2 | 3 ... 

• m = -l -{£-!) -(£-2) ...+(£-2) +(£-!) +£ forgiven £. 
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Notwithstanding their algebraic similarity (apparent from Table [I]), the 
integer and half-integer functions have different domains: the integer func- 
tions are defined on a Euclidean sphere, whereas the half-integer functions 
are not. One way of of interpreting this difference is that the angle has 
a range of O<0<47r for the half-integer functions compared with O<0<27r 
for the integer functions. This difference in domain is related to the essential 
difference between orbital angular momentum and spin angular momentum. 
The interpretation of the half-odd-integer functions is discussed in §[T2] on 
page 



12. 



Archival presentations of the associated Legendre functions [|], p. 332] in- 
volve generating functions and general formulae which in principle allow ex- 
plicit expressions and numerical values to be obtained for any required values 
of i and m . However some of these formulae become undefined when I and 
m are not integers, because they involve differentiation (with respect to 6) 
£ or m times; e.g. Rodrigues' Formula || §8.6.18 p. 334]. 

Thus to facilitate a clear exposition, we derive these functions as solutions 
of the appropriate differential equations. 

2 The 6 and (p Differential Equation 

The differential equation for the eigenfunctions, Y(6, <fi) and eigenvalues, A, 
of the square of the total angular momentum, L 2 is: 



.2 



1 ^4>1+^ctW#=^(m) 



l sin 00 V 36) sin 2 

It is well-known that the independent variables, and <fi , are separable, and 
hence Y(9, 0) can be written as a product: 

F(0,0) = 6(0)x$(0) (3) 

This leads to the separated, ordinary differential equations: 



sm0^ ( sin ^) +Asin 2 0-fi|e(0) =0 (4) 



and 

r d 2 



+ B $(0) = (5) 
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where B is the separation constant arising from the separation of 6 from 0. 
The solutions of are obvious (and well-known) to have the form: 

$(0) = exp(im0) (6) 

in which m is a constant to be determined. Substitution of (|6|) into (|5|) 
yields: 

(-m 2 + B) $(0) =0 or since $(0) ^ m 2 = B (7) 

This relationship between the exponent m and the separation constant B 
is necessary, but in itself it does not specify the solutions any further. It is 
noteworthy that B has the same value for the two different solutions of (|5|) 
having m values equal in magnitude but opposite in sign; e.g. for m=+2 
and m=— 2, B=4. 

In the case of orbital angular momentum one proceeds by noting that the 
angle takes any value within a complete circle (0 < < 2ir), and hence the 
appropriate condition is that the function $(0) shall have the same value 
when transits a complete circle [§, pp. 208-209], [§, p. 38]: 

$(0 + 2n) = $(0) exp(im[0 + 2tt}) = exp(im0) exp(im27r) = 1 (8) 

This requirement is necessary for $(0) to be a proper (i.e. single-valued) 
function of the points that comprise the surface of a sphere; it is met as long 
as m is any integer. 

We relax this single- valuedness condition on $(0) by leaving m undefined 
at this stage of the analysis, except that m should be real in order for $(0) 
to be periodic (and hence non-singular). Relaxing the traditional single- 
valuedness condition on $(0) , leads, however, to a different, non-classical 
interpretation of the angle (see § |372] on page p~2| ) . 



2.1 Solution of the 9 Equation 

Having solved the differential equation and determined that the separation 
constant B has the value B=m 2 , the 9 differential equation (|4]) (after dividing 
by sin 2 6) becomes: 



1 d ( . 

sin 6- 



sin#d# V d# 



A - 



2 

m 
shTfl 



= (9) 
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Following the usual derivation, we transform (|9]) to a new independent 
variable x defined by: 



x = cos 9 



sin 9 = y \ — x 2 



(10) 



and rename the dependent variable P(x); i.e. 0(9) = P(x). The range of 9, 
< 9 < 7i becomes +1 > x > — 1. 

With this transformation the differential equation (^) becomes: 



x 



d 2 P 

dx 2 



dP 
2x — + 

dx 



A 



m 



1-x 2 



P = 



(11) 



The last term, m 2 / (1 — a; 2 ), must be removed in order to develop the solution 
as a polynomial in x. This is achieved (as is well known) by writing the 
solution P(x) as the product of a known factor (1 — x 2 ) a and a to-be- 
determined factor, P'(x), that is anticipated to be a polynomial in x: 



P(x) = (l-x 2 ) a P'(x) (12) 
After substitution of the product form (pT2|) and its derivatives, the dif- 



ferential equation for P(x) (11) becomes an equivalent equation for the 
dependent variable P'(x): 



d 2 P' 

dx" 2 



'l-x< 



dP' ( 
- (2 + Aa)x— +[A-2a + 
dx \ 



Ax a — m 
1 — x 2 



P' = (13) 



The denominator (1 — x 2 ) will cancel out if we choose a as follows: 

^2 



4a z 



\m\ 



m 



a 



a 



and we choose the positive square- root: 

+ 



m\ 
~2~ 



(14) 



(15) 



because choosing a = — \m\/2 would make the solution (|1^) infinite at the 
ends of the range: x = ±1. 

With this choice for a, the differential equation ( |T3"D for P'(x) becomes: 



d 2 P' 
dx 2 



x 



2(\m\ + l)x 



dP' 

dx 
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[A- \m\(\m\ + 1)]P' = (16) 



We proceed to substitute a power series for P'( follows: 

i=0 



(17) 



and equating the coefficient of every power of x to zero (since the powers 
of x are linearly independent of each other) produces a set of homogeneous 
linear equations that determine the initial index k, the coefficients aj, and 
the eigenvalues of the separation constant A. 

Bypassing some details ||, it turns out that choosing k=0 produces all 
possible solutions, and thus we obtain the equation: 



i=0 



J2 (en [A - (|m| + i)((|m| + + a i+2 [(i + 2)(i + 1)] 



x 







18) 



Equating each power of x to zero produces a set of linear equations that 
determine the coefficients of the power series {cij : i = 0, 1, 2 . . .} which are 
concisely written as the matrix equation: 



/ T 00 
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V : J 




V ! / 



(19) 



which is summarized by the equation: 

T.a = 

where the elements of the square matrix T are given by: 

Tu = A — (\m\ + i)(\m\ + 1 + 1) T ii+2 = (i + l)(i + 2) 

2.2 The Eigenvalues of A 



(20) 



(21) 



Since the matrix T is triangular, its determinant is equal to the product 
of its diagonal elements, and hence the eigenvalues of A (denoted by Ai - 
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defined as being those values of A that make the determinant of T zero) 
are obtained by equating any one of the diagonal elements to zero: 

Ai = (|m| + i){\m\ + % + 1) for i = 0, 1, 2, . . . (22) 

This expression for the eigenvalues of the separation constant A shows that 
for a given value of |m| (for the two values of m: m = ±|m|) there is a set 
of values of A that increase quadratically with increasing values of i: 

A i = {\m\+i){\m\+i + l) =i 2 +i(2\m\ + l) + \m\{\m\ + l) £ = 0,1,2... (23) 

We see that the smallest value of Ai is \m\ (\m\ + 1) and that there is no upper 
limit on i nor upon the value of A4. The index i is the (mathematically) 
natural "quantum number" for designating an eigenvalue Af, the index i is 
the degree of the polynomial in x, P'(x), in equation (|T7|). 

However, it has become customary in physics to designate an alternative 
quantum number, £, defined by: 

£=\m\+i Ai = (|m|+i)(|m|+i + l)=^ + l) (24) 

£ = \m\, (\m\ + 1), (\m\ +2) ... 

Thus we see that £ can take any positive value beginning at |m| for a given 
value of the quantum number m in exp(im(j)). 

Furthermore it has also become customary in physics, to reverse the prece- 
dence of the relationship between £ and \m\ given in (J23J) by regarding £ 
as primary, with m as the secondary quantum number: 

£ = 0,1,2... A = £{£+!) 0<|m|<£ -£ < m < +£ (25) 

This physical viewpoint arises from atomic spectroscopy wherein the energy 
of an atomic state depends upon £, but all 2£+l states differing only in m 
are degenerate in the absence of an external magnetic field 0, p. 4]. 

For a given value of £ the minimum and maximum values of m correspond 
to m min = —£ and m max = +£. These 2£+l values of m correspond to 
polynomial degrees, i, differing by 1: 



% 





1 2 


2 


1 


m 


-£ 


-£ + 1 -£ + 2 . 


. £-2 


£-1 £ 
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and hence these 2£ + 1 values of m span an interval equal to 2£ that is 
symmetrical about zero. This interval is necessarily a non-negative integer, 
n , and hence: 

2£ = n a non— negative integer (26) 

and hence the allowed values of i are: 

I = \ (27) 

where n is any non-negative integer. When n is even, I is itself an integer, 
but when n is odd, i is a half of an odd-integer. This is how both integer 
and half-odd-integer values of I and m arise in the solution of the associated 
Legendre differential equation. This argument is essentially identical with 
that used in the abstract theory of angular momentum 0, pp.311] to deduce 
that both integer, and half-odd-integer, values of t and m are allowed. 

2.3 The Eigenfunct ions {a,i : i = 0, 1, 2 . . .} 

Choosing Tu = in ( pTf ) (i.e. Ai = [\m\ +i][\m\ + will lead to a solution 



in which the only non-zero coefficients are and all lower coefficients of 
the same parity (i.e. odd or even). The ratios of these non- zero coefficients 
are given by the recursion relation: 

[A-(\m\ + k)(\m\ + k + l))}a k + (k + l)(k + 2)a k+2 = (28) 
(k + l)(k 4- 2) 

" t = - (i-A)(2H +i + i + l) a ^ *»*=(i-2),(i-4)...{lar0} 

The recursion will terminate at k — 1 if z is odd, and at k — if z is 
even. Since k < i the factors in the denominator of (p8[) are always both 
positive, and hence consecutive terms of the power series in x alternate in 
sign. This recursion relation (p8| ) is written as beginning with the highest- 
index coefficient a^, from which the lower-index coefficients are calculated. 
Alternatively, the recursion can be written as beginning with ao or a±: 

(i - k)(2\m\ + i + k + 1) „ , . . 

= " (tll(t + 2) at *»* = 0orl...(i-2) (29) 
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Table 1: Legendre Functions Pj, m '(x) (£ = \m\ +i) 



\m\ 


factor 





1 


i 


= 2 


i = 


3 


i 


= 4 


i = 5 







1 




1 


x 


1 


-3a; 2 


3a;- 


-5x 3 


3 


-30a; 2 +35a; 4 


15a;— 


70x 3 +63x 5 


1 

2 


(1-x 2 ) 


i 

4 


1 


x 


1 


-4x 2 


3x- 


-6x 3 


3- 


-36x 2 +48x 4 


15x— 


80x 3 +80x 5 


1 


(1-a; 2 ) 


1 

2 


1 


x 


1 


-5a; 2 


3x- 


-7x 3 


3- 


-42a; 2 +63a; 4 


15a;— 


90a; 3 +99x 5 


3 
2 


(1-a; 2 ) 


4 


1 


x 


1 


-6x 2 


3x- 


-8x 3 


3- 


-48x 2 +80x 4 


15x— 


100x 3 +120x 5 


2 


(1-a; 2 ) 




1 


x 


1 


-7x 2 


3x- 


-9a; 3 


3 


-54a; 2 +99a; 4 


15a;— 


110a; 3 +143a; 5 


5 
2 


(1-x 2 ) 


5 

4 


1 


x 


1 


-8x 2 


3x- 


-10x 3 


3- 


-60a; 2 +120a; 4 


15a;— 


120a; 3 +168x 5 


3 


(1-x 2 ) 


3 
2 


1 


X 


1 


-9x 2 


3x- 


-llx 3 


3- 


-66a; 2 +143a; 4 


15a;— 


130a; 3 +195a; 5 


7 
2 


(1-x 2 ) 


7 
4 


1 


X 


1 


-10x 2 


3x- 


-12x 3 


3- 


-72x 2 +168x 4 


15a;— 


140x 3 +224x 5 


4 


(1-x 2 ) 


2 


1 


X 


1 


-llx 2 


3x- 


-13x 3 


3 


-78a; 2 +195x 4 


15a;— 


150a; 3 +255a; 5 


9 
2 


(1-x 2 ) 


9 
4 


1 


X 


1 


-12x 2 


3x- 


-14a; 3 


3- 


-84a; 2 +224a; 4 


15a;— 


160a; 3 +288a; 5 


5 


(1-x 2 ) 


■ r > 

2 


1 


X 


1 


-13a; 2 


3x- 


-15a; 3 


3- 


-90a; 2 +255a; 4 


15a;— 


170a; 3 +323a; 5 


11 
2 


(1-x 2 ) 


11 
4 


1 


X 


1 


-Ux 2 


3x- 


-16a; 3 


3- 


-96x 2 +288x 4 


15x— 


180x 3 +360x 5 



This recursion relation was used to generate Table the coefficient of x° or 
x 1 was chosen to be positive and all of the coefficients were re-normalized to 
make them the smallest set of integers for a given polynomial degree i. 

Table [1] displays a striking similarity between the well-known associated 
Legendre functions for integer values of i and \m\, and the newly discovered 
functions for half-odd-integer values. One difference between them, is that 
the factor (1 — x 2 )~ = sin' m ' 6 is an integer power of sin# when m is an 
integer, but that it has a factor of Vsin^ when m is a half-odd-integer; 
in the latter case the gradient, dPj m '(^)/d6 l is infinite at the limits: 0=0, 
9=7i, whereas in the integer case these gradients are finite. The polynomial, 
P'^ m \cos0), involves only integer powers of x = cos6 in both cases. 
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Table 2: Normalization Integrals Njj of Legendre Functions P' e m '(x) 



\m\ 


z 


1 


2 


3 


4 


5 





2 


2 


8 


8 


128 


128 




T 


3 


5 


7 


"sT 


TT 


i 








9tt 




257T 


2 


2 


8 


2 


32 


~2 




1 


4 


1 


32 


32 


768 


768 




3 


15 


21 


45 


"55" 


"9T 


3 


3tt 




15-7T 


3tt 


35-7T 


75tt 


2 


8 


16 


32 


16 


8 


32 


z 


16 


16 


64 


192 


6144 


2048 




15 


105 


45 


385 


455 


315 


5 


5tt 


5?r 


7tt 


35tt 


1357T 


945?r 


2 


16 


128 


16 


256 


32 


512 


3 


32 


32 


512 


512 


4096 


20480 




35 


315 


385 


1365 


315 


3927 


7 


357T 


7tt 


1057T 


27?r 


2079?r 


385vr 


2 


128 


256 


256 


256 


512 


256 


4 


256 


256 


1024 


1024 


16384 


81920 




315 


3465 


819 


3465 


1309 


19019 


9 


637T 


21vr 


99vr 


693tt 


10017T 


1287vr 


2 


256 


1024 


256 


8192 


256 


1024 


5 


512 


512 


4096 


4096 


32768 


32768 




693 


9009 


3465 


17017 


2717 


9009 


11 


2317T 


33tt 


3003vr 


1437T 


3861tt 


8775tt 


2 


1024 


2048 


8192 


2048 


1024 


8192 



2.4 Normalization of the Spherical Harmonics 

The normalization constant, iV 2 , of Y(9,<p) (equation |3|) is defined by: 



2tt 



N' 2 =l m<f))\ z d<f) I Q l ; n{ (6) ~ sm6d8 (30) 
J o Jo 1 1 

The (f) factor of iV 2 is always 2tt because |$(0)| 2 =exp(— imip) exp(im0) = l. 
However, since the domain of the functions for the half-odd-integer values 
of m, is O<0<4tt, the integration over involved in ([30]) should also, in 
principle, be over this range O<0<4tt. This would only make the <fi factor of 
iV 2 equal to 4tt rather than 2tt , and so for the sake of consistency we choose 
to have the same range of integration for both the integral and half-integral 
values of m. 



2 



10 



The 9 factor of iV 2 (Ng) is given by: 



P7T 




2 r +1 




2 




e^(9j 


sin9d9 = J 




dx 


Jo 









One difference between the case of integer values of £ , and that of half-odd- 
integer values, is that in the latter case Ng has a factor of tt, whereas in 
the integer case it is a rational fraction. This is illustrated in Table |2|, whose 
values correspond to normalizing the functions shown in Table [TJ 



3 Discussion 



The fermion quasi-spherical Harmonics were discovered by Dr. Ian Schlifer 
during a summer collaboration with Geoffrey Hunter several years ago ||. 
More recent collaborative work involving the other authors of this article led 
to the recognition of the potential utility of these functions for modeling the 
magnetic field of fermions such as the electron |], §]. 

The simplicity of the functions and their similarity to the well-known 
spherical harmonics for integer I and \m\ (Table p, suggested that they 
might have been discovered many years ago. However, a quite extensive 
literature search, including treatises on mathematical functions |J and on 
theoretical physics [fTPI , did not reveal any previous presentation of them. 



They have existed in principle as special cases of the hypergeometric 
function |, p.332 & pp.561-562], but: 

• the polynomial nature of their factors, P^'' m '(cos#), and 

• their potential application as eigenfunctions of the spin 
angular momentum of fermion particles, 

have not previously been recognized. 



3.1 Mathematical Aspects 

It is beyond the scope of this article to attempt a comprehensive considera- 
tion of these fermion spherical harmonics comparable with that in a standard 
treatise [|, pp. 332-341]. Here we simply point out some of the salient mathe- 
matical aspects to be considered in adapting the general theory of Legendre 
functions to the fermion quasi-spherical harmonics: 
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The formulae which are apparently undefined when £ and m are not 
integers (because they involve differentiation with respect to 9, £ or 
m times), notably Rodrigues' Formula §8.6.18 p. 334]: 

P o (r) 1 d"(x 2 -lf 

and the formula involving differentiation w.r.t. m ||, §8.6.6 p. 334]: 

dHp°(V) 



may be applicable by the theory of semi-differentiation jTT], p. 115 & 
p.307]. 

The recurrence relation |5|, §8.5 p. 334]: 

(£_| m | + l)pH = ( 2 £ + l)xPl ml -(e+\m\)P l ™l 



is applicable to the case of half-odd-integer values of £ and |m|; it is notewor- 
thy that all the coefficients in this recurrence relation: (£— \m\ + 1), (2£+ 1), 
(£ + \tti\), are integers even when £ and m are half-odd-integers. 

This recurrence relation produces a normalization and phase in which the 
highest power of x in each polynomial is positive, and in which the polynomial 
coefficients are generally fractions. A different normalization and phase was 
chosen to construct Table [I] in order to display the regularity of the series of 
functions. 



3.2 Interpretation 

The mathematical way of defining the half-odd-integer, quasi-spherical har- 
monics to be proper (i.e. single- valued) functions of the angle , is to define 
its range to be < <p < An, since: 

exp(im[0 + 47r]) = + exp(im0) 

when m is half of an odd integer. This is concordant with the well-known 
47T symmetry of fermion wavefunctions; i.e. the angle must transit two 
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complete circles for the wavefunction to return to its original value [O, p. 21 



& p.138], H p.141]. 

However, with the range of redefined in this way, <fi can no longer 
be regarded as one of the coordinate angles of the points on the surface of a 
sphere. Rather it is an angular coordinate of the points on a double-sphere, 
for which the points on the outer surface of the sphere are different from 



the corresponding points on the inner surface [16], p. 419]. The geometry 
and topology of the double sphere can be modeled by the Dirac belt trick 
fj~6| , p.417ff]. Such a closed-on- itself surface with 4tt periodicity is the well- 
known Mobius band |TJ], pp. 141-143], p6| , p. 418]. Thus the fermion spherical 
harmonics are not true spherical harmonics, but rather quasi- spherical har- 
monics. 

An alternative way of dealing with these functions is to define their do- 
main as the spherical range 0<4><2tt, but then to recognize that they are 
double-valued functions of 0. This is the approach of Bethe's theory of 
double-groups as discussed by Altmann |L3|, chapter 13]. It has also been 
based upon the theory of non-simply connected spaces [|14], §3.2 pp. 35-38]. 

In group-theoretical terms the integer i and m functions are representa- 
tions of £0(3) while the half-integers functions are representations of SU(2) 



14], pp.35-41], 0, pp.123-130], 0, pp.140-142]. 
Regardless of whether one takes: 

• the single-valued function on a non-spherical domain approach, or 

• the double-valued function on a spherical domain approach, 

the angle <fi does not have a classical, physical interpretation. Notwith- 
standing the long history of spin angular momentum J7L its physical nature 
remains something that is not easily interpreted in terms of rotational motion 
in classical space-time. 
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